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1 . INTRODUCTION 


^ This  report  validates  and  documents  the  operation  of 
the  Multi-Doppler  Quadrature  Correlator  (MDQC) . A previous 
technical  note*  described  the  analytic  technique  utilized  in 
the  development  of  the  MDQC.  ^ 

Section  2 of  this  report  describes  the  method  used 
to  validate  the  results  obtained  from  the  MDQC.  Included  in 
this  section  are  the  expected  theoretical  values  and  a summary 
of  the  actual  results.  The  computer  program  used  in  the 
validation  is  found  in  Appendix  A. 

Appendices  B,  C,  D,  and  E contain  descriptions  and 
program  listings  of  the  special  TIMFAX  Black  Boxes  required  in 
order  to  use  the  Multi- Doppler  Quadrature  Correlator. 


TRACOR  Document  No.  SD/69-008-U,  "A  Technique  for  a Multi- 
Doppler  Quadrature  Correlation  Processing  Scheme,"  5 August  1969. 
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2.  VALIDATION  OF  THE  MULTI -DOPPLER 
QUADRATURE  CORRELATOR  (MDQC) 

The  operation  of  the  Multi-Doppler  Quadrature  Correlator 
was  validated  by  using  the  MDQC  routine  to  correlate  a linear  FM 
slide.  Because  of  the  range-Doppler  ambiguity  of  an  LFM  pulse, 
correlation  peaks  are  obtained  in  all  channels.  The  relative 
delay,  At, between  adjacent  channels  depends  on  the  bandwidth,  B, 
and  pulse  length,  T,  of  the  signal,  and  is  given  by 

At  - Af  T/B,  (1) 

where  Af  is  the  channel  spacing.  For  the  LFM  up-slide,  which  was 
used  in  the  validation  program,  the  peaks  in  the  up-Doppler  channels 
will  occur  later  than  the  ones  in  the  zero-Doppler  channel.  The 
peaks  in  the  down-Doppler  channels  will  occur  earlier. 

Since  neither  the  signal  nor  the  reference  have  been 
time-compressed,  there  is  no  dilation  loss  involved  here.  The 
amplitude  of  the  correlation  peaks,  relative  to  the  zero  Doppler 
channel,  is  determined  by  the  overlap  loss,  and  is  given  by 

(2) 

The  validation  program  consists  of  the  TIMFAX  topology 
shown  in  Figure  1.  The  Black  Box,  PGENl,  generates  an  LFM  slide 
from  100  to  200  Hz,  sampled  at  1000  Hz.  The  pulse  is  1 second 
long.  CORGNl  and  C0RGN2  generate,  respectively,  sine  and  cosine 
representations  of  the  same  LFM  signal.  The  references  from  the 
two  uses  of  C0RGN  are  each  broken  into  7 pieces  by  the  two  uses 
of  MS CAT.  MS CAT  also  clips  the  references,  and  reformats  them 
into  a form  suitable  for  MDQC.  The  7 output  channels  from  MDQC 
are  processed  through  7 uses  of  LOCAL  and  7 uses  of  LOCPR  to 
extract  and  print  the  maximum  value  in  each  channel  and  the  time 
at  which  it  occurred. 

MDQC1  uses  an  8 point  Fast  Fourier  Transform  to  produce 
7 Doppler  channels.  The  spacing  between  the  channels  is  thus  7/8 


At  _ -I 

Af 

T 

B 
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Hz.  The  effect  of  low-pass  filtering  prior  to  comb  filtering 
which  is  inherent  in  the  design  of  MDQC  is  given  in  Equation  12 
of  Reference  1 as 


I 


u)  T 
sin  7F" 

«>oT 


(3) 


where  oj0  is  the  Doppler  offset  (in  radians/second),  T is  the  pulse 
length  in  seconds  and  N is  the  number  of  pieces  in  the  reference. 
Since  the  channel  spacing  is  7/8  Hz, 

ou0  = 2tt  -g  n,  (4) 


where  n is  the  channel  number,  measured  from  the  zero-Doppler 
channel.  The  filter  loss  on  the  nth  channel  is  thus 


I 


n 


sin 


rrn 

F" 


rrn 

F" 


(Note  that  this  factor  can  usually  be  ignored.  However,  the  test 
case  violates  the  restriction  given  by  Equation  1 of  Reference  1, 
resulting  in  higher  than  normal  filter  loss.) 

The  results  of  the  validation  program  are  predicted  to 
be  as  follows: 

(a)  The  amplitude  of  the  zero  Doppler  channel  for  a 
clipped  reference  will  be  equal  to  the  average 
absolute  value  of  the  input  signal.  Since  PGENl 
produces  a sinusoid  with  a peak  amplitude  of  1, 
the  peak  in  the  zero-Doppler  channel  will  be  2/tt 
» .6366. 

(b)  The  amplitude  of  the  nth  channel  will  be  modified 
by  Equations  2 and  5 and  will  be 


3 


(c)  The  time  of  the  peak  in  the  nth  channel  will  be 


cn  " 1 + mi  ’ (7) 

where  n > 0 represents  the  nth  up-Doppler  channel 
and  n < 0 represents  the  nth  down-Doppler  channel. 

The  results  predicted  by  Equations  6 and  7 are  tabulated 
in  TABLE  I,  together  with  the  actual  results  obtained  from  exer- 
cising the  MDQC  routine.  The  agreement  is  excellent. 

TABLE  I 

COMPARISON  OF  THEORETICAL  VALUE  WITH 
COMPUTED  RESULTS  FROM  MDQC 


Channel 

Measured 

Amplitude 

Theoretical 

Amplitude 

Measured 

Time 

Theoretical 

Time 

-3 

.4868 

.4861 

.973 

.9737 

-2 

.5619 

.5631 

.983 

.9825 

-1 

.6171 

.6149 

.991 

.9913 

■ 0 

.6359 

.6366 

.999 

1.0000 

1 

.6139 

.6199 

1.008 

1.0088 

' 

2 

.5690 

.5631 

1.016 

1.0175 

3 

.4957 

.4861 

1.025 

1.0263 

; 

im 

1 

« 

4 
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APPENDIX  A 

VALIDATION  COMPUTER  PROGRAM 
FOR  MDQC 
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APPENDIX  A 

VALIDATION  COMPUTER  PROGRAM  FOR  MDQC 
TIMFAX  TOPOLOGY 
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APPENDIX  B 

MDQC  BLACK  BOX  DESCRIPTION 


PROGRAM  NAME: 
NATURE  OF  WORK: 
PROGRAMMER'S  NAME 
DATE: 

COMPUTER  LANGUAGE 


MDQC 

Black  Box 
David  Doan 
October  25,  1969 
FORTRAN  V 


PROBLEM  DEFINITION 

MDQC  is  a Multi-Doppler  Quadrature  Correlator.  The 
Cooley-Tuckey  Fast  Fourier  Transform  is  used  as  a comb  filter  to 
produce  the  effect  of  a number  of  frequency  shifted  channels. 

The  output  of  each  filter  is  processed  through  an  ideal  envelope 
detector. 


OPERATING  INSTRUCTIONS 

1 . Topology  Card: 

The  topology  card  requires  four  input  channels  as 

follows : 

MDQCi  MSCATj  ,MSCATk,FFTIN<t  ,S0URCEm 

The  first  two  fields  must  be  produced  by  MSCAT;  the  first  one  is 
the  real  part  of  the  reference,  and  the  second,  the  imaginary 
part.  The  third  field,  which  must  be  produced  by  FFTIN , contains 
the  cosine  reference  required  by  the  Fast  Fourier  Transform. 

The  time  function,  indicated  above  by  SOURCEm,  is  the  input 
signal  channel. 

2 . Data  Input: 

All  P-array  entries  are  fixed  point  unless  noted. 

P(l)  - Number  of  outputs  (teeth)  from  comb  filter 
P(2)  = Number  of  output  channels,  s P(l) 


1 


B-l 


P(3)  = Size  of  Fast  Fourier  Transform, 

! P(l) 

^ P(l)  of  MS CAT 

= P(l)  of  FFTIN 

. 

= a power  of  2 

P(4)  ■ Output  scale  factor,  SV  ■ 1.0  (floating  point) 

P(5)  ■ Estimated  length  of  correlator  reference  * 
actual  length 

P(6)  - Estimated  input  record  length  ? actual  record 
length 

P(7)  * Number  of  pieces  into  which  reference  is  broken. 

SV  provides  automatic  transfer  of  this  value  from 
MS CAT.  See  Note  below.  Must  equal  P(l)  of  MSCaT. 

P(8)  = 0 for  fixed  point  time  function  input 

* 1 for  floating  point  time  function  input. 

SV  provides  automatic  transfer  of  this  value 
from  MSCAT.  See  Note  below.  Must  equal  P(3) 
of  MSCAT. 

NOTE:  If  the  fields  from  MSCAT  are  input  directly  to 

MDQC,  value  for  P(7)  and  P(8)  need  not  be 
supplied.  If  the  fields  are  being  read  from  the 
drum  or  tape,  P(7)  and  P(8)  must  be  supplied. 

See  Section  8,  Program  Conditions  and  limitations  for 

a further  discussion  of  the  P-array  values  required. 

3 . Input  Fields: 

The  input  fields  must  be  produced  by  the  boxes  shown 
in  (1)  above.  The  topology  for  generation  of  the  complex  refer- 
ences should  be  one  or  the  other  of  the  following  to  provide  for 
Doppler  shift  which  increases  with  channel  number. 


TF 
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I 


f 

I 


If  the  references  are  directly  generated,  as  for 
instance  with  CORGN,  then  the  topology  should  be 


05  C0RGN1 
MSCATl 
CORGN 2 
MSCAT2 


CORGN 1 


sine  wave 


C0RGN2 


cosine  wave 


03 

MDQC1  MSCATl, MSCAT2 ,FFTIN1 , . . . 

If  the  references  are  to  be  obtained  by  quadrature 
heterodyning,  the  topology  should  be 


13 

- - - 

PGEN1 

cosine 

PGEN2 

sine 

TMUL1 

PGENl , SOURCEn 

TMUL2 

PGEN2 ,S0URCEn 

FLTPL1 

FSETl ,TMUL1 

FLTPL2 

FSETl ,TMUL2 

RFEXT1 

FLTPLlB 

RFEXT2 

FLTPL2B 

05 

MSCATl 

RFEXTl 

MSCAT2 

RFEXT2 

03 

- - - 

MDQC1 

MSCATl ,MSCAT2 ,FFTIN1 , . . 

Output  Fields: 

Channel  A is  used  as  a scratch  field.  It  should  not 
be  referenced  by  any  other  Box.  It  can  be  overlayed  by  fields 
not  used  in  the  Section  3 containing  MDQC. 

5 • Input  Time  Function: 

The  input  time  function  may  be  either  fixed  or  floating 
point.  The  mode  requested  is  determined  by  P(8). 


i 
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6.  Output  Time  Functions: 

The  output  time  functions  are  floating  point.  P(2) 
channels  are  produced.  If  P(l)  > P(2),  the  channels  are  block- 
multiplexed  in  the  following  manner: 

Channels  1 through  P(2)  - 1 are  multiplexed  by  a factor 
of 


M 


1 


Int 


PCI)  + P(2)  - 1 

PIT) 


Channel  P(2)  is  multiplexed  by  a factor  of 

M2  = P(l)  - (P (2 ) - 1)  • M1  . 

Note  that  some  combinations  of  P(l)  and  P(2)  will  cause  M2  to 
be  zero.  This  situation  should  be  avoided,  since  extra  storage 
is  reserved  which  will  not  be  used. 

Block  multiplexing  means  that  a full  record  (equal  in 
length  to  the  input  record)  is  output  for  each  multiplexed 
channel,  followed  by  a full  record  of  the  next  multiplexed 
channel.  This  is  in  contrast  to  the  sample  by  sample  multiplexing 
used  by  MPX  and  DEMPX . The  block-multiplexed  output  must  be 
demultiplexed  by  BDMPX. 

If  the  output  is  to  be  held  in  floating  point,  P(2) 
must  be  less  than  or  equal  to  17. 

The  first  time  function  output  is  on  Channel  B.  The 
Doppler  shift  associated  with  each  output  from  the  filter  is 

. ) 

TOT  r pml 

f = I n - Int 

P(3)T  L 1 -1 

where  P (1)  is  the  P(l)  input  to  MSCAT,  f is , the  Doppler  shift 
in  Hertz,  T is  the  pulse  length  in  seconds,  and  n is  the  "tooth" 
number  prior  to  multiplexing.  P(3)  in  this  expression  will  be  a 
power  of  2.  If  P(l)  is  even,  there  will  be  one  more  up-Doppler 
channels  than  there  are  down-Doppler  channels. 
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The  format  of  the  output  time  function  is  as  follows: 


Channel  B 


Channel  C 


Last 

Channel 
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will  be  produced,  indicating  that  P(3)  has  to  be  set  to  the  next 
higher  power  of  2. 

The  following  messages  indicate  fatal  errors.  The  run 
is  aborted  if  any  of  them  appear. 

MDQCn  UNBALANCED  REFERENCES 

NO.  OF  PIECES  j,k 

T-F  MODE  l,m 

The  number  of  pieces  or  the  time  function  mode  was  not  the  same 
for  both  references. 

MDQCn  UNBALANCED  REFERENCES 

SHORT  REF  LENGTH  ■ j ,k 
NO.  LONG  REFERENCES  l,m 

The  length  of  the  two  references  was  not  the  same. 

MDQCn  OVERFLOW  IN  p*  STORAGE,  PIECE  j 

LAST  REQUEST  = k,  1 available 

Probably  caused  by  too  small  a value  for  P(6),  the  estimated 
record  length. 

MDQCn  P(3)  = IS  ILLEGAL 

is  produced  if  P(3)  < P(l) 

MDQCn  INCREASE  P(5)  BY  AT  LEAST  

P(5)  is  less  than  the  total  length  of  the  reference  supplied. 

MDQCn  PARAMETER  ERROR 

EXPECTED  PIECES  = , REAL  = , IMAG  - 

EXPECTED  MODE  = , REAL  = , IMAG  * 

P(7)  is  not  equal  to  P(l)  of  MSCAT  or 
P(8)  is  not  equal  to  P(3)  of  MSCAT. 


J 


W"«wu 


frwxmul 

AHHHnCVTM 
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MDQCn  WRONG  FFTIN  FIELD 
EXPECTED,  SUPPLIED 

P(3)  does  not  agree  with  P(l)  of  FFTIN 

THE  FIELD  FROM  OR  

IS  AN  ILLEGAL  INPUT  TO  MDQCn  UNLESS  P(7)  AND 
P(8)  ARE  MANUALLY  SET. 

If  the  references  are  recorded  on  tape  or  drum,  the 
parameters  P(7)  and  P(8)  must  be  supplied  in  the  data  input. 
If  they  are  not  supplied,  the  first  two  input  fields  must  be 
supplied  directly  from  MS CAT. 


The  following  diagnostics  were  included  for  check-out 
purposes  and  should  not  occur: 

MDQCn  CORR  RECORD  TOO  LONG,  PIECE  

MDQCn  REFERENCE  ERROR,  I = 

NRR  = , REALF  = 

NIR  = , IMAGF  = 

ERROR,  REFERENCE  LENGTH  SUPPLIED  TO  MDQCn  IS  SMALLER 
THAN  THE  LENGTH  OF  THE  REFERENCE  SUPPLIED 

8.  Program  Limitations  or  Conditions: 

MDQC  must  be  in  Section  3. 

The  parameters  for  MDQC  and  MS CAT  are  as  follows: 
(P-array  values  are  for  MDQC,  unless  indicated  by  an  over-score, 
Over-scored  values  are  for  MSCAT.) 

a.  Determine  the  maximum  Doppler  shift,  fmax,  from 
other  data  known  about  the  problem  (relative 
velocities,  range  of  interest,  etc.).  This  value 
may  not  exceed 
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ym 

fmax_^  = for  1 dB  dilation  loss,  or  (1) 

fmax_3  = — for  3 dB  dilation  loss,  (2) 

where  fo  is  the  center  frequency  of  the  signal  in 
the  water,  B is  the  bandwidth,  and  T is  the  pulse 
length. 

b.  The  number  of  pieces  is  limited  by 

P(l)  i 4 fmax  T.  (3) 

c.  The  frequency  resolution  is  given  by 


P7TT 

Af 

P(3)  T . (4) 

To  minimize  splitting  loss, 

Af  * "ZT  * or  (5) 

P (3)  a 2 PHT  . (6) 

d.  P(3)  must  be  a power  of  2. 

e.  Juggle  P(l)  and  P(3)  to  satisfy  Equations  (3)  and 
(6),  and  to  provide  a convenient  value  of  Af. 

f.  The  total  number  of  outputs  teeth  is  then 


P(1)  = _2_fmax_T^.(32  + x 

TOT) 


(7) 
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The  input  references  must  be  constructed  in  accordance 
with  the  description  in  Section  3,  above.  Otherwise,  the  sign 
of  the  Doppler  shift  may  be  changed. 

WARNING:  If  P(l)  and  P(2)  are  such  that  the  time  function  outputs 
are  multiplexed,  the  outputs  must  be  demultiplexed  by  using 
BDMPX  before  writing  the  output  to  tape  or  drum.  This  requirement 
exists  because  the  output  record  length  from  a tape  - or  drum  - 
read  Box  is  not  necessarily  rhe  same  as  the  record  length  at  the 
input  to  the  tape  - or  drum  - write  Box.  Under  these  conditions, 
the  subsequent  demultiplex  operation  will  be  hopelessly  confused. 

9.  Storage: 

2771g  or  152910  words 
10.  Subroutines  Referenced: 

C0RR2,  C0RR3,  FIXTST,  FLTTST,  TW0P0W,  BBNAME , PVAL, 
CALSAL,  CTFFT,  EREXIT,  NPRT$ , NI01$ , NI02$,  SQRT,  NERR3$ , CBNAME. 

PROCESSING  METHOD 

MDQC  is  a digital  simulation  of  the  type  of  correlator 
in  which  a comb  filter  follows  the  multiplier.  The  general 
arrangement  of  such  a correlator  is  shown  in  Figure  1. 


Input 


time 

comb 

compressor 

— -0— 

filter 

memory 

bank 

Doppler 

Outputs 


reference 

memory 


Figure  1 
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There  are  two  problems  in  directly  simulating  the 
simple  method  shown  in  Figure  1.  The  sampling  rate  at  the  input 
to  the  comb  filter,  which  is  equal  to  the  input  sampling  rate 
times  the  time  compression  factor,  is  extremely  high.  This  makes 
simulation  of  the  filter  very  time  consuming.  The  method  shown 
is  unable  to  distinguish  between  positive  and  negative  Doppler, 
since  the  pass  band  of  an  ordinary  filter  depends  only  on  the 
absolute  value  of  the  frequency. 

These  problems  are  solved  in  the  following  way: 

(1)  The  reference  is  replaced  by  a quadrature  (complex) 
reference . 

(2)  The  multiplier  is  followed  by  an  "integrate-and- 
dump"  low  pass  filter. 

(3)  The  comb  filter  is  simulated  by  using  the  Fast 
Fourier  Transform  (FFT) . 

This  arrangement  is  shown  in  Figure  2. 


Input 


Doppler 

Outputs 


Figure  2 
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The  integrate-arid-dump  low  pass  filter  greatly  reduces 
the  sampling  rate  at  the  FFT  ("comb  filter"),  thus  reducing  the 
running  time.  The  use  of  a quadrature  reference  in  combination 
with  the  Fast  Fourier  Transform  provices  a method  for  retaining 
Doppler  sense. 

Mathematically,  the  Multi- Doppler  Quadrature  Correlator 
consists  of  the  evaluation  of  the  equations  given  below. 

The  output  of  the  low  pass  filter  is  given  by 


N 

K 


K(n-KL) 

XT 


Xk+m  rm 


m = 


Kn 

N 


(1) 


where  N = the  number  of  pieces  into  which  the  reference  is  broken 
K = the  number  of  samples  in  the  complete  reference 
* the  signal  input  (real) 
amd  rm  = the  complex  reference. 


The  output  of  the  comb  filter  is  given  by 

2rr-tn 


p(k,i)  = 


V1 


v X I»(k)e'3 


N 


n = 0 


(2) 


which  is  iust  the  discrete  Fourier  Transform  of  I (k) . N is  a 
J n o 

power  of  2,  greater  than  or  equal  to  N. 

The  parameter,  k,  is  the  sample  number  of  the  output 
time  function.  The  parameter,-  i,  is  the  Doppler  parameter.  It 
is  related  to  the  Doppler  shift  at  the  center  frequency  by 
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where  T is  the  length  of  the  pulse,  and  may  be  either  positive 
or  negative. 

The  absolute  value  indicated  in  Equation  (2)  performs 
the  function  of  an  ideal  envelope  detector.  No  further  rectifica 
tion  or  averaging  is  required. 

For  a more  complete  mathematical  treatment  of  this 
process,  see  D.  B.  Doan,  "A  Technique  for  a Multi-Doppler, 
Quadrature  Correlation  Processing  Scheme,"  TRACOR  Document  No. 
SD/69-008-U,  August  1969. 
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APPENDIX  C 


MDQC  PROGRAM  LISTING 


subroutine:  mdqco  (n,p> 

OIMENSIOM  p(  I > .'REALFU  >V»  M AGF  < 1 ) , A I (I  )~»  COREF ( I , ,PC(  1 i , JOuTTI  lT» 

I N 0 UT  (JJ  ,A  J ( 1 ) , N AM  t ( 2 J NAME ( 2 ) , KNAMEC2) 

implicit  Integer  (a~-Y> 

real  B,_  A l,  A U i P C ,__C 0 RE  F j__  .Z 

parameter  P T Gj  = 8 

DEFINE.  Q(  I )=P(  I +P T Q_) 

C*  T h E V A K I A 0 L E (J  l II  IS  THAT  PORTION  0F~  THE  P-aKRAY  USED 

C FOR  INTERNAL  SCRATCH  STORAGE 

dittn a at  mn rsx  it o 

DEFINE  AAJ(  1 )»AJ(  l ) 

DEF  I NE  PpC  < 1 ) sPC ( I ) 

EQUIVALENCE  < NRR  ,NAME(  II  I * (NXR»NAmE(2)  ) 

E X T L K N AL  CO HR~2  ,"‘c 0 R R 3 
N=  1 3 

p(‘Mseool(  i .0)  « std  value  for  scale  factok 

PI7)s-_l  6 5V  FOR  AUTOMATIC 

P ( a'lV-'l  "(3  sv  AKA  MET"  e H ■ T R AN  ST  E ft 

RETURN 

entry  MDQC  i TiTi'npfVi nptT mnpt'» l as ,"i outT,  iYue'.  i ou t TkoUT t7a~»b7 

CALL  FIXTST  I P ( 1 ) > P ( 3 ) ) 

Call  FLTTST  ( P ( *» ) »P  ( H j ) 

CALL  FIXTST  <P(5)jP<a>) 

N C T s f >*  0 P 0 irs  ( P ( 3 ) i 5 720  ) (3  IS  P ( 3 ) a POivER  OF  2 
IFINCT  .EQ.  P ( 3 ) } GO  JO  5 51  YES 
_C  ALL  "CBN  A HE  '(N  A ME NO PRINT  MSg'aNDSET  f6"~N~EXT ’’hTsHER 

PRINT  IOQ'4,  NAME,  NCT  (3 POWER  OF  2 

P (j  f=NCT 

5 CALL  CBNAHE  J 1 , JNAME ) 

CALL  CBN A ME  ( 2 ,’KNAmE)  " 

IF  ( AMD(  JNAME  ( 1 ) -SHMSCAT  , -63  ) _•  EQ  * 0 .AND. 

i and'T'Kname ( i )”  " s”h hs c at,”  -63  "7  ‘.TqV TTgo  To"  s' 

IF  (P(7)  « G T » 0 .AND.  P(6)  ,6E.  0)  GO  TO  7 U J,,MP  IF  NOT  AUTO 

call  BONaME  (NAMfr) 

PRINT  1012, JNAMEjKNAME, NAME 
CALL  ERE X I T 

6 P ( 7 I =PVAL  ( 1 . 1)  ffl  NO,  PIECES  FROM  REAL  REF  GENERATOk 
p(ft)=PVALU  .3)  M T-F  MODE  FROM  R E A L ' RE  F " GENE  K Ai~OK 

I F ( P ( 7 ) ,N£.  P V A L ( 2 i 1 ) ) GO  TO  700  y ERROR  IF  RfcAL  AND  IhAo 

fF  (~P  (V)  . N £ • P V AL  ( 2 » ^ ) ) 60  TO  700  fl  PTRAMS  NqT  Tht  SAME 

7 I F ( P I 7 ) ,GT.  P ( 3 ) ) GO  TO.  720  0 FFT  AREA  TOO  SMaLL 

- ,NpfV3  

INPTsl 

j-NP  T Bp  raT" 

LAG*P<5)  * SS 


IttunHiiH 
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toutf»i 

g(2)e|P(|)  ♦ P(2)  - U/P<2)  lo>  OUTPUT  MPX  FACTOk 
0 o i»  ( 2 ) 

Q(j)sli(2)  • P ( 2 ) IS  PIECES  AHlCH  CAN  BE  STORED  IN  Aj 
N S » 2 • P < 7 ) w total  pieces 

LC»P(5)/P(7)  0 PREDICTED  SHORT  PEP 

LC  C = N OD ( P ( 5)  , > ( / ) ) Id"  NO  . LONG  REFS  ' 

<uh)-o  .s  counter  for  pieces  in  ai  

NKh*0  13  COUNTER  FOR  AI  AREA  USED 


NIRsO  ® COUNTER  FOR  AI __AREA  AVAILABLE  _ 

do  ill— i“mTns  r«i  Coop  to  allocate"  sYo'h'age 

l F { I .LE.  0(3))  GO  TO  9 » JUMP  IF  IN  AJ 

N R R «N R R +"  P ( 6 ) B P ( 6 ) I S EST I M ATED  I NPUT  RtCORo”  LENG  TH  " 

I F l NR  R . G T . N I R ) GO  TO  IS  0 JUMP  I F_.T  00  MUC_H 

oT4')»Q(V)'M  m count  pieces"  in  ai  ‘ 

9 IF  1 ANO(  1,1)  . N E . 0)  GO  TO  10 

N I R a N 1 R + LC  Q RELEASE  AFTER  EACH  PAIR 

I F ( ( 1 + 1 ) /2  . LE  . LCC  ) N I R 3 N I R + 1 _ fl_  I NCR  FOR  LOiiG  R_tF 

10  'CONTINUE"'  s'  END'LL L 0 Cat  ION  LOOP 

IS  NCT  = NS-0(3_)-U(S)  _W__P_IECES  LEFT  T0_  STORE 

TfTn'cT'VlTT'o')  'N'cTsO  ' S PREVENT  ACCIDENTS 

I S I ZE  = P ( 7 ) * 1 08  + P ( 3 ) * 2 + NCT*P(6>  P SCRATCH  FlELD 

(T(  fa  f=  I s I i E IS  SIZE  SCRATCH  F'lELO  n 

CALL  B B N A ME  (NAME) 

pr i n't"  i'Qo» ," name , "p <y+p' TQ")7“PTN ♦ p t gl  V YctT  nsV'Tsue"' 

«ct)=Q(N)  + 0(3)  fj  total  pieces  in  a i_  and  aj_ 

fCMJTTlT)VpT2T'3  NOV  QU PUT  "CHANNELS 

I0uTT(2)3l  p MAY  BE  UNBALANCED 

KOuTTsl 

A s 0 

RETURN 

entry  mdqo  i_P) 
q u r-I'VfMTf  ia'lTz  a’tTon  ' f l a'g~ 

RETURN 

ENTRY  MDRC3  (P,KKLG,REALF,IMAGF»C0REFiIM,INPT,aJ,LaG,PC,JI1, 

1 10 U T T ,_A  J »_NOUT  ) ___ 

lag  WorTeT  sample's 

IF  (JMC  ,GT.  0)  GO  TO  13 
JMCVO 

J IsD 

‘(TO  T 0 I Hb 

lb  IF(Qll)  . NE • 1)  GO  TO  30  W SKIP  INITIALIZATION 

— rr»-RrAXFTirrvNE'."p'rr)  i . e orrery  act *r par", a keteft  £■*  ko  k 

IF  ( K E A t.  K ( S ) .NE.  P<8)>  GO  TO  760  0 PARAMETER  ERROR 
i ft  i"m  A GFT T)"'Vn  e'.">T7TT  "g  0"'TO"  7 6"0"'ra  " p" A R A Ml  E Y L R"'e"r  « 0 « 

IF ( I M A G F ( h ) .NF,  P ( 8 ) ) GO  TO  760  M PARAMETER  ERROR 

ITT R E A L r f2~)  TSI ET.  I'M  a GF ( 2 H G 0 TO  7 cTl  W ER R 0 R Ip  KtA L a'NO  I MAG 

I F < R E A L F ( 3 ) ,NE.  IMaGF(3))  GO  TO  701  0 REF  L c.  N G T H S ARE  NOT  EmUAL 

NCTsCOREF  ( 2 T"S"ChECK'TENGTH'  Of'-  FFT'  REF  

I F ( P ( 3 ) .NE.  NCT)  GO  TO  750  0 ABORT  IF  DOESN’T  MATCH 
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NS  a"2'*  P<7)  13  TOTAL  PIECES  IN  tjOTH'REFS 

NCT=0  W P ARRAY  POINTER 

N l<  K b7~  U _K  E F E R FN  C E~  ARNA  Y~PO  I NTER 

Z = BOOL(P(4)  )/FL0aT(REALF(6)  ) 

00  20  'I  »TVNS' " 

PC < NCT+  I ) =UOOL < RE ALF < 2 ) ) G P(l)  FOR  COrR,  KLF  tENGTh 

IFt  ( T+  ir/2  . LEV  RE ALF  ( 3 ) ) PC  ( NCT+  1)  *800L  ( BOOL  ( PCY^CTVI  ) ) ♦ jT'ieT'LO  N'G 

PC(NCT+2)«FL0AT(KEALF(NRR+1 ) )*Z  f,l  P ( 2 ) FOR  CoKR,  SCALt 

P CTN Cr^l  sRO 0LTPT3 T)  13  "~PT3  ) ~F OR~COrRV  "T-F  11OI1E 

CALL  CALSAV  ( 1 ,C0RR2,PC(NCT+l ) ) » INITIALIZE  Co«R 

TFI  ANOCr,!  ) ‘VEQ  ♦“"QT'URRaNRR  “V  REA  L F(N  R RfVTi,'"I  N CTT"  RE  F POINTER'"’ 

20  N C T = N C T ♦ 5 4 0 INCR  P ARRAY  POINTER 
q ( 1 > =o  X>  CLEAR  'I  N I T I AL  I 2 AT  I ON">XAG 
30  LC=LAG  - P ( 7 ) *RE ALF  ( 2 ) - REALF<3)  a BASE  LAG  V/.LUE 
I F1XX~.  L T~ SV)~ G O T O 7 3 0 
NC  T a l (3  OUTPUT  STORAGE  SWITCH 
NCP*!'  W'P'  ARRAY ‘ POINTER 
N R K a 7 M REAL  REF  POINTER 

N I R'a  7~F3  I'M  A G~R  E F~P  0 I NT  E R 

N S = P ( 7 ) ra  NO,  PIECES  IN  EACH  REFERENCE 
N 0 C T a r M -'A  J-S 7 OR  aGE"PO  I NTER 

N I C T = 1 W A I STORAGE  pointer 

NSCfaPl  7 8 ~P  C~  ST  OR  A G E P 0 I NTER 

00  100  1 = 1, NS  13  begin  correlation  LOOP 
LCCaLOl 

L C = L C + REAL F (2)  a add  rep  length  to  lag 

1 F ( I • L E , RE  ALF  (3  ) ) LC  = LC+’  1 0 INCR  IP  LONG  Rtr 

JNaJMC  ♦ LC  W .TOTAL  . SAMPLES.  TO  _CORR 

IF  I NC  T « G T • N ( 3 ) J GO  TO'  40 

CALL  CALSAV  <11,C0RR3,PC(NCP),KFLG,KEALF<NkK),jN,1,A1,LC,JMC,1, 

s XJ  ( U O CTT7N  O : 

_ NOCT=NOCT  + JMC  « INCR  A J STORAGE  POINTER 

GO” TO  60  

SO  IF ( NC T ,gT,  Q<4))  GO  TO  50 

Ca  l c "c  a ls’a  v"(TiVC'o'r  r'3  Vpc'(  ncpT,  k Flg  , real  fin  rkTi  i n ,T',  a'iTlcT  j'MCVlV 
5 AI(NICT),NC) 

MICTaTJTCT  ♦ Jifc  ra  IN'CR  A I STORAGE  POINTER 
IF  ( N I C T ,GT.  LCC)  GO  TO  702  « A_I  OVERFLOW 

GO  TO  6 0 * 

50  CONTINUE 

CA  LL"  C A L S A V-fH  VC0RR3 VP C I NC PT.  K F XG 7* E ALF ( N R R 7 y ‘AT, XC  V jT-TCV l '• 

S PCINSCTJ.NC) 

N S C T a N S C T ♦ JtlC  W INCR  PC  STURAGE  POINTER 

IFIMSCT  ,GT.  0(5)  “ 2 * P 1 3 ) ♦ 1 > GO  TO  703  a P_C__OVEKFLO». 

"VO  "I  F t NC  » G T •'  J M C ) ~6 0 ~ T 6 706  ~7s ~C  ORR~~OU TPU  T RECORD  TO 0 L ON G 

NR R aN.RR.  + REA LF  (NR R I_.«.  _I  NC _K  _ .R.E_A.L__RE_F._Pp_l.N_T t_K_  

"TfTnRK  ,G  T.  REAL  nil'VTlio  TO  740 
NCTaNCT  ♦ I 0 INCR  STORAGE  SWITCH 
NCPatJCP  ♦ b 4 « I NCR  P ARRAY  POINTER 
IP iNCT  » 6T » 0(3))  GO  TO  70 


■■MHMMMMHMi MMMW> 


C-3 
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call  calsav  J iT,'cork'3»pc  < ncfT.k fl  g“,'im a g F < n i k ) , jn,  i . a j ,lc , jmc,T*‘ 

$ A J ( NOCT ),NC) 

NOCtakOtT  ♦"jHC"  U"I  N CR  A J STORAGE  P O'fN TER  

GO  TO  9 U 

~T fi  T T~f  N "CT  TG  T~.  QT‘0-)  G 0 T 0 B 0 

call  calsav  ( i i ,corr3,pc<ncp) ,kflg, imagfini*) , |N, i , a j ,lC,jmc,  i » 

5 Al  (NICT)  ,NC>  *■ : 

NICT  = MCT  + JMC  w I NCR  Al  STORAGE  POINTER 

If  (NICT  , G T • ' LCC)  GO  TO  704  A T OVERFLOW"  ~ 

GO  TO  90 
8 0 C O n TT  N U E 

CALL  CALSAV  ( l 1 , C0RR3 , PC < NCP  ) , KFLG , I M AGF ( N I R ) , | N , i , A I ,LC, JMC  , 1 . 
S'  pc  (NSCT  ) . r J C ) 

NSCT  = NSCT  ♦ JNC  INCR  PC  STORAGE  POINTER 

IF  (NSC  T .GT.  0(5)  - 2*P(3)  ♦ I)  GO  TO  705  i*'PC  OVE"KF"lOv. 

90  1F<NC  • GT » JMC)  GO  TO  706  li  OUTPUT  RECORD  FROM  CORK  TOO  LONG 
W IK  = NTR  + I"Ma  G nriTRl't'l  I'NCR  I'M  A G REF  PQ  I NT  £.  R 

I F ( N I R .GT.  IMAGF(l)  ♦ 2)  GO  TO  7HQ 

NCT  = NCT  + 1 0)  INCR  STORAGE  SWITCH 

l QO  NCPoNCP  + S4  w INCR  p ARRAY  POINTER 

END  OF  CORRELATION  LOOP 

NSs2*P(7)-l  ra  TOTAL  PIECES  MINUS  1 
— N SC  Taw  ( 5 ) - 2*  PT3l~ B~U  A S E “OF  F F T "~W  0 R K TNG  S T OTYacYe 

NF1=NSCT*1  (3  LIMITS  FOR  CLEARING  FpT  STORAGE  J 

N F 2 - 0 ( 5 ) 

N I C T = P ( I ) - I « MO.  TEETH  IN  FILTER  MINUS  1 

HTR  =PTn7‘2~r"o F’F SETT  FOR- N EG  aTTVE  F REQUEn'C  IES 

do  iso  1=1, jmc  ® sample  loop 

DO  1 05  " J = NF  1 , NF  2 

105  PCIJJ-0  U clear  fft  WORKING  storage 

NCT-NSCT  Q INITIaLI/E  fft  input  pointer  " 

DO  130  J = 0 , NS  .3  P ILCE  LOOP 

CT  ♦ 1 « INCR  FFT  INPUT  POINTER 
.GE.  W ( 3 ) ) GO  T_0_  l 10 

P C ( N C T > = A A~J~(T  Vj  MC  ) ©GET  SAM  P L E FROM  a J 
GO  JO  _1  30 

iT6"r  fTj  "rer; " wrrs  o'-tirw 

PC(NCT)  = aAM I + (J-Q(3)  ) » JMC  ) S GET  SAMPLE  FROM  ul 
GO  TO  130 

20  P C ( N C T ) « PP  C U ♦ ( J - Q ( 4 ))  * J M C ♦ P ( 7 )_♦_!  0 a _)  0 G E T S AM  p_Lj  E_J_  R 0_M__P  C 
I"3  0-~C0"ntT  NO  L ‘ 

CALL  CTFFT  ( PC  < NSCT+ 1 ) , CORE F (4)  , P ( 3 ) » 2 > 

”0  0 l"4  0 ZTs  07 NT CT"TY"  0 uTPUT-"ST0RA  G 6r"L  OOP 

NOCT  = J - NIK  0 INDEX  OFFSET  FOR  NEGATIVE  F K E G u L_N  C_|_E_S 

IF(NOCT  .LT.  OJNOCTsNOCT  ♦ P ( 3 ) 

1 A A J < 1 .7J?jJLlC-La-S£RJ_iP_P£Lii5Pl*2*_Np_CT^lJ_*_*2__t__PP_C_(jJiS_Qjjt2*_N0tJ_+2J_»_*j£_>. 
ENO  OUTPUT  LOOP 

1 45  NCTalOUTTU  ) __  

0 oTs  0 f«  2,NC  T fj  SET  N0~»  6 U T P U T~S  A MPL  E S 
ISO  N 0 U T ( I - 1 ) ° J M 


/ 
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' 


i/t- 


\ 


I 


NOUT ( NC T ) = < P ( 1 ) - ( P ( 2 ) - 1 > »Q < 2 ) ) • JMC 

RETURN 

/On"  N C T * P V A L_(  2 , 1 )” 

NCPsPVAL (2,3) 

pk  i n t ' i uoOV  p (7  iV'NCrr  p < 8 )7"NCP  

GO  TO  7 V 9 

7 oT~ PR  IN  T 10  0 1 , REALFT2)  » I~M  A GF  ( 2),  REALM  3),  IMAgF(3) 


GO  TO  799 

7 02~  N 0 C T = ‘IH REAL  

GO  TO  710 

' — 7 0 2TN  0 C T = 4 H R E A L 

GO  TO  711 

7Xm~~ti  a CT  = 4 H 1 T17VG 

GO  TO  7 1 n 

705"N0CT  = 4HIM AG 

GO  TO  711 

~J-7a6"PR  I NT  1003V 1 

GO  TO  799 

7T0  NR R = NIC  T - I 

N I K = L C C - 1 

-NjCT  = 2HAT  

GO  TO  715 

7Yl " N R R a|Ni S C T '-‘"I 

N I R = CJ  ( 5 ) - 2 • P ( 3 ) 

N fCT=2HPC 

715  PRINT  1002,  NlCT,  NOCT,  I,  NRR,  NlR 

GO  TO  799 

720.  PRINT  l JOG  , P j(_3  P.C  7J 

7 3 U N C T = b 5 - LC 

PRINT  1 U U 6 . tiCJ 

GO  TO  799 

7 4 0 PRINT  1 Qo  7 , I , _NRR  RE_ALF_U  Lj_.NJ_R.t_J  «AGf  t_U_ 

GO  TO’  799 

750  PRINT  10n9,  P ( 3 ) , N C T ; 

GO”  TO  7 9 9’ 

760  PRINT  1010,  P(7),  REA  L F ( 4^  , I M A G F ( 4 ) , _^P  ( »__n  t A LF  ( & ) t 

"s  l MAGF ( 5 ) 

• GO  TO  79v  

79  9 Call  wBNaiTE  (name) 

PRINT  IQ  1 1 . NAME : 

CALL  ekexit 

I LOG  F okmaT  ( l X , 2_x_ , _'_lj N_B _A_L  A_N C_£ P_ _ RJE F_E R_E N_C E _S ' /_  * N0_* Q g P I E££ S * J_b /_*  _ _T "_F 

SMOOE  * bX , 2 I b ) 

1001  FORMAT  ( t X , 1 2X  . »UN  0 A iL  AN_C  £0  RE  F E RE INCE.S!  /_!_ 

SO.  "'LONG”  REFS*  ,4X,2l5)” 

1032  FORMAT  ( 1 X , 1 2 X , * OVF.KKLOv.  IN  « A 3 , « ST  OR  AGE  , * Ab  , , P 1 ECL  ♦ , 1 5/ « LAST  RE 

WES  T"=~*TsT*  » *Tb  , ♦“ TV  A IL  aTTL E * ) 

1303  FO R M AT  ( 1 X , 1 2 X , * C 0 R R RECORD  T 00  L On  G , P I EC E * j l b JL 

" 1 004  'FORM aT(  IX ,2A6,’P ( 3 ) "C HA N G E 0 TO  *14) 

1035  FORMAT  ( l X , 1 2 X , ♦ P < 3 ) ■ • , I 4 t * IS  ILLEGAL. FOR  * l 4 ,*PI  ECES*  ) 


I 


i 


C-5 


- " — 


L 


1 006 
1 007 

1 008 


1009 

1010 

ToTT 
.101  2 
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FORMAT  ( I X , 12n  , « INCREASE  P(S)  OY  AT  LEAST  Mb) 

JFO  KM  AT  (J  X ,J2  X • ’REFER LNC  E ERROR,  l = M5/»  NKK»»|S,»,  RE ALF ( I ) ■ M b/ « 

i N1R«* IS,* , I MA6F ( 1 )»♦ 1 5 ) 

FORMAT  ( 1X.2A6, ’PIECES  STORED  IN  Aj»’l5,’,  AI»iJb,’,  SCRATCH.  M S 
»:7  T 0 T AL>T£  C l S = •":&,  V,"  SCR  AT  CM  'STOR  A CiE"  M sf~~ 

FORMAT  < t X , l 2 X , ’ V R 0 N G FFTIN  F I ELD_,_»  I & , • E X PEC T f 0 , • I & , * SUPPLIED*) 

format  (YxVi^xV* parameter'  ekrorV/  » e x p e cYe  b">Vt  c e s* *T& e a l ■ *T&Y 

i«  , IMA G=«IS/»  EXPECTED  HOPE  »M5,i,  RE  ALa  * I S , « , IP.Ab»*lb) 

F 0 R M AT  ( 1 M + , 2 A 6 ) 

FO  R M AT.  J_  t T HE FJ_£  L D_  F R 0 M„!A6jAU  .*._M__»  A.*jJU_il 

• IS  AN  ILLEGAL  INPUT  To  »A6,A1, 

! * _ 0 N L E_S_S_ _P _<_7 )_ _ AN _D_ _ P_( 8>_ _ A_ R_£_ _ _M_  A_N  U A_U kIL  _S  ET  •__*__  1 _ 

END 


I 

I 

[ 

! 

I 
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APPENDIX  D 

MSCAT  BLACK  BOX  DESCRIPTION 


PROGRAM  NAME: 
NATURE  OF  WORK: 
PROGRAMMER'S  NAME 
DATE: 

COMPUTER  LANGUAGE 


MSCAT 
BLACK  BOX 
David  Doan 
August  25,  1969 
FORTRAN  V 


PROBLEM  DEFINITION 

MSCAT  produces  a set  of  partial  correlator  references 
from  a single  input  field  for  use  with  the  Box,  MDQC.  The 
technique  used  by  the  Black  Box,  SCAT,  is  used  to  provide  very 
efficient  clipped-reference  correlation. 


OPERATING  INSTRUCTIONS 

1.  Topology  Card: 

The  topology  card  for  MSCAT  requires  one  floating 
point  field  such  as  that  provided  by  CORGN  or  RFEXT. 

2 . Data  Input: 

All  P-array  entries  are  fixed  point. 

P(l)  = number  of  pieces  into  which  the  reference  is 

to  be  broken.  See  Section  8 of  the  MDQC  write-up 

P(2)  = length  of  output  field.  Initially,  P(2)  should 
be  equal  to  the  length  of  the  input  field.  It 
may  be  possible  to  reduce  this  value  to  conserve 
storage.  See  Sections  5 and  6. 

P(3)  = O(SV)  for  fixed  point  time  function  input  to  MDQC 
= 1 for  floating  point  input  to  MDQC. 

3.  Input  Fields: 

The  input  field  must  contain  a floating  point  correlator 
reference  in  the  format  produced  by  Boxes  such  as  RFEXT  or  CORGN. 


D-l 
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4.  Output  Field: 

The  output  field  contains  the  count  word,  and  five 
constants,  followed  by  P(l)  subfields,  each  having  the  format 
of  the  field  generated  by  SCAT.  The  field  is  shown  below: 

word  1 
word  2 
word  3 
word  4 
word  5 
word  6 
word  7 
word  8 

word  N-^ 
word  N-^ 
word  N-^ 

word 

etc. 

5 . Printer 

One  lim 
The  format  of  the 

MSCATx 

where  x is  the  use  number,  N is  the  number  of  samples  in  this 
partial  reference,  K,  L,  and  S are  parameters  indicating  the 
operating  of  SCAT  (See  SCAT  write-up  and  SCAT/CORR  process 
description),  and  y is  the  number  of  the  partial  reference. 
Following  the  above  is  a line  with  the  following  format: 

MSCATx  FIELD  STORAGE  AVAILABLE  = , USED  = 

The  storage  available  is  P(2)  + 5.  The  storage  used  is  the  total 
number  of  words  in  the  output  field  including  the  count  word. 


+ 7 
+ 8 
+ 9 


+ 1?2  + 9 


Output: 


N 

length  of  short  reference 
number  of  long  references 

P<1) 

P (3) 

total  reference  length 
N1 

constants  and  machine  instructions 
as  generated  by  SCAT  for  piece  1 


N2 

constants  and  machine  instructions 
as  generated  by  SCAT  for  piece  2 


2 is  printed  for  each  piece  of  the  reference, 
line  is 


N = K = , L = , S = , PIECE  y 
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6 . Program  Limitations  and  Conditions: 

a.  Each  piece  must  contain  at  least  2 samples, 

b.  An  error  message  is  produced  if  the  output  field 
is  too  small  and  the  run  is  aborted. 

7 . Storage  Requirements: 

447g  or  2951q  words . 

8.  Subroutines  Referenced: 

SCAT3,  FIXTST,  CALSAV,  BBNAME,  EREXIT,  NPRT$ , NI01$ , 
NI02$ , NERR3$ . 

PROCESSING  METHOD 

MSCAT  breaks  the  reference  supplied  in  the  input  field 
into  P(l)  pieces.  If  the  input  field  contains  N samples,  then 
there  are 

n = N mod  P(l) 

long  pieces  and  N-n  short  pieces. 

The  first  n pieces  have  a length 

11  = Int  (N/P (1) ) +1, 

and  the  remaining  N-n  pieces  have  a length 

12  = Int  (N/P (1) ) . 

Each  piece  is  converted  to  a series  of  add-subtract 
commands  by  SCAT3  (part  of  the  Black  Box,  SCAT) . The  outputs 
from  SCAT3  are  stored  in  the  output  field. 


D-3 
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APPENDIX  E 


MS CAT  PROGRAM  LISTINGS 


II 


SUliKCUT  l N E MSCATO  |N»P) 

DIMENSION  Hill,"  REFdV,  0UT<1),  NAmE(2)" 

EQUIVALENCE  tNS.NAMEU)),  ( OCT  , NAME ( 2 ) ) 

I'm  Pl  fC  I T I N T E G E R ( A - Z ) 

RtAL  KEF , SAVE 
E X T E K N A L SC  A T 3 " 

N 3 6 

P ( 3 ) * 0 ~ a ~ F I X ED  " T » F ~ IS  "STANDARD 
RETURN 

TivTr-y— f-nsrcXr  1 (PTfNpr,  iITpTTkTnpt , l a a , lOWmifUE,  i uu  fr , KOUT  T » A • 0 ) 
I n p y = i 
1 0 U T E = 1 ~ 

I S I ZE  = P ( ? ) ♦ S 

call  fixtSt  “(>(  i iVpUVj 
RETURN 

E >.“7RY~M  S ~CAT~2  fP) 


RETURN 

E r.T R Y MS C A T 3 ( P , XF LG  » RE F , Ou T ) 

NCT=REr(|)  a LENGTH  OF  REF  SUPPLIED 
OUTl  2)  = NC  T / P (~f  )""~M  ” L £ N G T H 0~F  SHOP  T ~R  E F 
Ojj  T ( 3 ) = i-l_00j NCT  ,_P ( 1 ) ) g NO*  0 F LONG  REFS 
0 U T(4  1 = P ( 1 ) E1  NO,  OF  PI  EC  E’s 
OUT  ( S ) =P  ( 3 } fi  T-F  MODE 
OUT  { 6 ) «NC  T »'  TOT  At  "REFER  E"NCE  LEn  GTH 

P(S)«P(2I  a P(l)  For  SCAT,  length  of  ref 
p < b ) -.T"¥"PT2)  FOR  SCAT,  FLO  aYVng  "point""fulo^>  l"a  g 

p < A ) = P ( 3 ) Q P ( 3 ) FOR  SCAT  , T-F  FLAG 

NCT* 1 © INPUT  FIELD  P 0 INTER 
H f,  T = 7_  fs»_  OUTPUT  _F_lE_LD__  POINTER 

ns * >TiT  «'"n oT'of  pieces" 

DO  1C  1 = 1 »_N S _ 

OC  T=0UTT2  ) l’e  NG t"h""6 f'"S H 0 R T""R E F " 

1_F  ( I . L E_.  OUT  ( 3 ) ) 0 CTaQC T*J_  jrf_  ] N_C R F 0_R_  L ON G R E (■ 

SaVLbREF  (NCT  ) « SAVE  oTrIO  CONTENTS  OF  INPUT  fl£.LO 
REF«NCT)=OcT  a SET  REF  LENGTH  IN  FIELD  FOR  PAkTIALKEF 

call  calsav  i *» ,scat3 , p ( 4 1 Ykflg ,R E F("nC T ) ,OUTj MC f ) ") 

PRINT  1000  , I I?  ADD  PIECE  NO.  TO  SCAT  MSG 
Rtf  ( HC T ) «S  A V E ~ G>  R E S T 0 RE  I N P UT  F I E L D 
NCT  = NCT  ♦ OCT  lal  I NCR  INPUT  POINTER 


l 


I 


E-l 


BMP 


/ 
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TM 


"To  mct  = kct  + ouumcti'  ♦ l Q inch  output  poi'ntlr  ~ 

If ( NC  T - 6 .GT,  P ( 2 ) ) GO  TO  30  g OUTPUT  FIELD  T&o  SMALL,  crfkGft 
■ OUT  l I ) = M C T - 2 9 SET  LENGTH  Of  OUTPUTTiELD 

2 C N C T c P ( 2 ) +50  STORAGE  AVAILABLE 

MCT=MCT  - V 0 STORAGE  USED 

call  b b n a m e (name) 

prTnt  lUOl,  NAME",  NCT,  MCt 

I f ( 0 U T ( 1 ) ,LT,  0)  CALL  EWE  X I T 0 ABORT  RUN 

--Return-" 

3d  PRINT  1002 

OUT  (T)  * - T cf  SET  E R R 0 ft  FLAG 
GO  TO  20 

-IlTDI)  T OTTHa '!  ' < f hTTTS“XT*T"I  E"Cr*TTT) 

luUl  FORMAT  ( lX,2AiS,  'FIELD  STORAGE  A V A I L A OLE* * , l b , • . uS£D«*,Ib) 

10  02  F 0 ft  M A T " ( / * ’•••ERROR***'*  ) " 

iT  *.i  *"\ 


